Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm 
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We present the first results of the Fermilab Muon g—2 Experiment for the positive muon magnetic 
anomaly a, = (gą —2)/2. The anomaly is determined from the precision measurements of two angu- 
lar frequencies. Intensity variation of high-energy positrons from muon decays directly encodes the 
difference frequency wa between the spin-precession and cyclotron frequencies for polarized muons 
in a magnetic storage ring. The storage ring magnetic field is measured using nuclear magnetic reso- 
nance probes calibrated in terms of the equivalent proton spin precession frequency à; in a spherical 
water sample at 34.7?C. The ratio wa /@,, together with known fundamental constants, determines 
a4 (FNAL) = 116592 040(54) x 1071} (0.46 ppm). The result is 3.3 standard deviations greater than 
the standard model prediction and is in excellent agreement with the previous Brookhaven National 
Laboratory (BNL) E821 measurement. After combination with previous measurements of both pt 
and L^, the new experimental average of a,(Exp) = 116592061(41) x 10^! (0.35 ppm) increases 
the tension between experiment and theory to 4.2 standard deviations. 


INTRODUCTION 


The magnetic moments of the electron and muon 


Be = ge (5) S8 where ge =2(1 + aj), 
2m, 


(£ = e, p) have played an important role in the develop- 
ment of the standard model (SM). One of the triumphs 
of the Dirac equation [1] was its prediction for the elec- 
tron that ge — 2. Motivated in part by anomalies in 
the hyperfine structure of hydrogen [2, 3], Schwinger [4] 
proposed an additional contribution to the electron mag- 
netic moment from a radiative correction, predicting the 
anomaly [5] ae = «/2-» ~ 0.00116 in agreement with 
experiment [6]. 


The first muon spin rotation experiment that observed 
parity violation in muon decay [7] determined that, to 
within 10%, g, = 2, which was subsequently measured 
with higher precision [8]. A more precise experiment [9] 
confirmed Schwinger’s prediction for the muon anomaly 
and thereby established for the first time the notion that 
a muon behaved like a heavy electron in a magnetic field. 
This evidence, combined with the discovery of the muon 
neutrino [10], pointed to the generational structure of the 
SM. 


The SM contributions to the muon anomaly, as illus- 
trated in Fig. 1, include electromagnetic, strong, and 
weak interactions that arise from virtual effects involv- 
ing photons, leptons, hadrons, and the W, Z, and 
Higgs bosons [11]. Recently, the international theory 
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FIG. 1. Feynman diagrams of representative SM contribu- 
tions to the muon anomaly. From left to right: first-order 
QED and weak processes, leading-order hadronic (H) vacuum 
polarization and hadronic light-by-light contributions. 


community published a comprehensive [12] SM predic- 
tion [13] for the muon anomaly, finding a,(SM) = 
116591 810(43) x 1071! (0.37 ppm). It is based on state- 
of-the-art evaluations of the contributions from quan- 
tum electrodynamics (QED) to tenth order [14, 15], 
hadronic vacuum polarization [16-22], hadronic light-by- 
light [11, 23-36], and electroweak processes [37—41]. 

The measurement of a,, has become increasingly pre- 
cise through a series of innovations employed by three 
experimental campaigns at CERN [42-44] and more re- 
cently at Brookhaven (BNL E821) [45]. The BNL net 
result, with its 0.54 ppm precision, is larger than a, (SM) 
by 3.7 standard deviations (c). While the electron mag- 
netic anomaly has been measured to fractions of a part 
per billion [46], the relative contribution of virtual heavy 
particles in many cases scales as (m,/me)? ~ 43,000. 
This is the case e.g. for the W and Z bosons of the SM 
and many hypothetical new particles, and it gives the 
muon anomaly a significant advantage when searching 
for effects of new heavy physics. Because the BNL re- 
sult hints at physics not included in the SM, Experiment 
E989 [47] at Fermilab was constructed to independently 
confirm or refute that finding. In this paper, we report 
our first result with a precision of 0.46 ppm. Separate pa- 
pers provide analysis details on the muon precession [48], 
the beam dynamics corrections [49], and the magnetic 
field [50] determination. 


EXPERIMENTAL METHOD 


The experiment follows the BNL concept [45] and uses 
the same 1.45 T superconducting storage ring (SR) mag- 
net [51], but it benefits from substantial improvements. 
These include a 2.5 times improved magnetic field intrin- 
sic uniformity, detailed beam storage simulations, and 
state-of-the-art tracking, calorimetry, and field metrology 
for the measurement of the beam properties, precession 
frequency, and magnetic field [47]. 

The Fermilab Muon Campus delivers 16 highly po- 
larized, 3.1GeV/c, ~120ns long positive muon beam 
bunches every 1.4s into the SR. A fast pulsed-kicker mag- 
net deflects the muon bunch into a 9-cm-diameter storage 
aperture, resulting in ~ 5000 stored muons per fill. The 


central orbit has a radius of Ro = 7.112m and the cy- 
clotron period is 149.2ns. Four sections of electrostatic 
quadrupole (ESQ) plates provide weak focusing for ver- 
tical confinement. 

The muon spins precess in the magnetic field at a rate 
greater than the cyclotron frequency. The anomalous 
precession frequency [52] in the presence of the electric 
E and magnetic B fields of the SR is 


B ql 
a, B — day (- i1 


1 B x E 
(o ps i) c 
For horizontally circulating muons in a vertical magnetic 
field, B .B = 0; this condition is approximately met in 
our SR. At the muon central momentum po, set such that 
yu = y (1 + 1/aj) ~ 29.3, the third term vanishes. 

In-vacuum straw tracker stations located at azimuthal 
angle @ = 180° and 270? with respect to the injec- 
tion point provide nondestructive, time-in-fill dependent 
beam profiles M(a,y,¢,t) by extrapolation of decay 
positron trajectories to their upstream radial tangency 
points within the storage aperture [53]. These profiles 
determine the betatron oscillation parameters necessary 
for beam dynamics corrections and the precession data 
fits discussed below. 

Twenty-four calorimeters [54-56], each containing a 
9 x 6 array of PbF2 crystals, detect the inward-spiraling 
decay positrons. When a signal in a silicon photomul- 
tiplier (SiPM) viewing any crystal exceeds ~ 50 MeV, 
the data-acquisition system stores the 54 waveforms from 
that calorimeter in a set time window around the event. 
Decay positron hit times and energies are derived from 
reconstruction of the waveforms. 

The magnetic field is measured using pulsed proton 
NMR, calibrated in terms of the equivalent precession 
frequency v, (1;.) of a proton shielded in a spherical sam- 
ple of water at a reference temperature T, = 34.7?C. 
'The magnetic field B is determined from the precession 
frequency and shielded proton magnetic moment, Hp (T) 
using hw, = 25, B. The muon anomaly can then be ob- 
tained from [57] 
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where our collaboration measures the two quantities to 
form the ratio 


(3) 


The Run-1 data, collected in 2018, are grouped into 
four subsets (a — d) that are distinguished by unique 
kicker and ESQ voltage combinations. The ratio Rj, can 


be conceptually written in terms of measured quantities 
and corrections as 


felock wa (1 + Ce + Cp + Cmi + Coa) 
Feats (wh (2, y, $) x M(z y, 9)) (1+ Be + By) 
(4) 

The numerator includes the master clock unblinding fac- 
tor felock, the measured precession frequency wp, and 
four beam-dynamics corrections, C;. We deconstruct 
Gy» (T) into the absolute NMR calibration procedure (in- 
dicated by fci») and the field maps, which are weighted 
by the detected positrons and the muon distribution av- 
eraged over several timescales ((w,,(x,y, 6) x M (v, y, 9))). 
'The result must be corrected for two fast magnetic tran- 
sients B; that are synchronized to the injection. 

Damage to two of the 32 ESQ high-voltage resistors 
was discovered after completion of Run-1. This led to 
slower-than-designed charging of one of the quadrupole 
sections, spoiling the symmetry of the electric field early 
in each fill. The impact of this is accounted for in the 
analysis presented. Brief summaries of the terms in Eq. 4 
follow. 


" 
RF 


ANOMALOUS PRECESSION FREQUENCY 


fuod: A single 10 MHz, GPS-disciplined master clock 
drives both the wa and n measurements. The clock has 
a one-week Allan deviation [58] of 1 ppt. Two frequen- 
cies derived from this clock provide the 61.74 MHz field 
reference and a blinded “(40 — €) MHz” used for the wa 
precession measurement. A blinding factor in the range 
+25 ppm was set and monitored by individuals external 
to our collaboration. felock is the unblinding conversion 
factor; its uncertainty is negligible. 

we: The signature of muon spin precession stems from 
parity violation in u* decay, which correlates the muon 
spin and the positron emission directions in the pt rest 
frame. When boosted to the lab frame, this correlation 
modulates the e* energy (E) spectrum at the relative 
precession frequency Wa between the muon spin and mo- 
mentum directions. The rate of detected positrons with 
E > Eg as a function of time t into the muon fill then 
varies as 


N(t) = Nony (t)e re 
x [1 + Ana(t) cos (Wat + Yo + n9 (t))), (5) 


where yr, is the time-dilated muon lifetime (~ 64.4 us), 
No is the normalization, A is the average weak-decay 
asymmetry, and qo is the ensemble average phase angle 
at injection. The latter three parameters all depend on 
Fin. The n; terms model effects from betatron oscilla- 
tions of the beam, and are not required in their absence. 
'This beam motion couples with detector acceptance to 
modulate the rate and the average energy, and hence the 


average asymmetry and phase, at specific frequencies. 
The coherent betatron oscillation (CBO) in the radial 
direction dominates the modulation. 

The CBO, aliased vertical width (VW), and vertical 
mean ((y)) frequencies are well measured, and the mj 
terms are well modeled and minimally correlated in fits 
for wa. 

An accurate fit to the data also requires accounting 
for the continuous loss of muons over a fill, also weakly 
coupled to wa. Coincident minimum-ionizing energies in 
three sequential calorimeters provide a signal to deter- 
mine the time dependence of muon losses. 

Two complementary reconstruction algorithms trans- 
form the digitized SiPM waveforms into positron ener- 
gies and arrival times. In the “local” approach, wave- 
forms are template-fit to identify all pulses in each crys- 
tal, which are then clustered based on a time window. 
In the “global” approach, waveforms in a 3 x 3 array 
of crystals centered on a local maximum in time and 
position are template-fit simultaneously. After subtrac- 
tion of the fit from the waveforms, that algorithm it- 
erates to test for any missed pulses from multiparticle 
pileup. To avoid biasing wa, we stabilize the calorime- 
ter energy measurement within a muon fill by correcting 
the energy reconstruction algorithm on the SiPM pixel 
recovery timescale (up to tens of nanoseconds) and the 
fill timescale (700 us) using a laser-based monitoring sys- 
tem [59]. The system also provides long-term (many- 
days) gain corrections. The two reconstructed positron 
samples are used in four independent extractions of wa in 
which each e* contribution to the time series is weighted 
by its energy-dependent asymmetry; this is the optimal 
approach [60]. Seven other determinations using addi- 
tional methods agree well [48]. Each time series is modi- 
fied to statistically correct for contributions of unresolved 
pileup clusters that result from multiple positrons prox- 
imate in space and time. The analyses employ one of 
three data-driven techniques to correct for pileup, which 
would otherwise bias wa. 

A x? minimization of the data model of Eq. 5 to the 
reconstructed time series determines the measured (m) 
quantity w™. The model fits the data well (see inset 
to Fig. 2), producing reduced x?s consistent with unity. 
Fourier transforms of the fit residuals show no unmodeled 
frequency components, see Fig. 2. Without the 7; terms 
and the muon loss function in the model, strong signals 
emerge in the residuals at expected frequencies. 

The dominant systematic uncertainties on Wa arise 
from uncertainties in the pileup and gain correction fac- 
tors, the modeling of the functional form of the CBO de- 
coherence, and in the wcgo(t) model. Scans varying the 
fit start and stop times and across individual calorimeter 
stations showed no significant variation in any of the four 
run groups [48]. 

The measured frequency w” requires four corrections, 
Ci, for interpretation as the anomalous precession fre- 
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FIG. 2. Fourier transform of the residuals from a time-series 
fit following Eq. 5 but neglecting betatron motion and muon 
loss (red dashed); and from the full fit (black). The peaks 
correspond to the neglected betatron frequencies and muon 
loss. Inset: Asymmetry weighted e* time spectrum (black) 
from the Run-1c run group fit with the full fit function (red) 
overlaid. 


quency wa of Eq. 2. The details are found in Ref. [49]. 

Ce: The electric-field correction C, from the last 
term in Eq. 1 depends on the distribution of equilib- 
rium radii rz, = x — Ro, which translates to the muon 
beam momentum distribution via Ap/po = z«(1—n)/ Ro, 
where n is the field index determined by the ESQ volt- 
age [49]. A Fourier analysis [49, 61] of the decoherence 
rate of the incoming bunched beam as measured by the 
calorimeters provides the momentum distribution and 
determines the mean equilibrium radius (xe) ~ 6mm 
and the width coz, ~ 9mm. The final correction factor is 
C, = 2n(1 — n)8?(z2)/ R$, where (x2) = 02, + (ae)?. 

Cy: A pitch correction C, is required to account for 
the vertical betatron oscillations that lead to a nonzero 
average value of the B - B term in Eq. 1. The expres- 
sion C, = n(A7)/AR$ determines the pitch correction 
factor [49, 62]. The acceptance-corrected vertical ampli- 
tude A, distribution in the above expression is measured 
by the trackers. 

Extensive simulations determined the uncertainties 
dC. and dC, arising from the geometry and alignment 
of the plates, as well as their voltage uncertainties and 
nonlinearities. The nonuniform kicker time profile ap- 
plied to the finite-length incoming muon bunch results in 
a correlation introducing the largest uncertainty on Ce. 

CQ: Any bias in the average phase of muons that 
are lost compared to those that remain stored creates a 
time dependence to the phase factor yp in Eq. 5. Beam- 
line simulations predict a phase-momentum correlation 
dyo/dp = (—10.0 + 1.6) mrad/(%Ap/po) and losses are 
known to be momentum dependent. We verified the cor- 
relation by fitting precession data from short runs in 
which the storage ring magnetic field, and thus the cen- 
tral stored momentum po, varied by +0.67% compared to 
its nominal setting. Next, we measured the relative rates 


of muon loss (ml) versus momentum in dedicated runs 
in which muon distributions were heavily biased toward 
high or low momenta using upstream collimators. Cou- 
pling the measured rate of muon loss in Run-1 to these 
two correlation factors determines the correction factor 
Cmi. 

Cpa: The phase term qo in Eq. 5 depends on the muon 
decay coordinate (x,y,$) and positron energy, but the 
precession frequency wa does not. If the stored muon 
average transverse distribution and the detector gains 
are stable throughout a fill, that average phase remains 
constant. The two damaged resistors in the ESQ sys- 
tem caused slow changes to the muon distribution during 
the first ~ 100 us of the measuring period. An exten- 
sive study of this effect involved: a) generation of phase, 
asymmetry, and acceptance maps for each calorimeter as 
a function of muon decay coordinate and positron en- 
ergy from simulations utilizing our GEANT-based model 
of the ring (gm2ringsim); b) extraction of the time de- 
pendence of the optical lattice around the ring from the 
COSY simulation package and gm2ringsim; c) folding the 
azimuthal beam distribution derived from tracker and 
optics simulations with the phase, asymmetry, and ac- 
ceptance maps to determine a net effective phase shift 
versus time-in-fill, yo(t); and d) application of this time- 
dependent phase shift to precession data fits to deter- 
mine the phase-acceptance (pa) correction Cpa. The use 
of multiple approaches confirmed the conclusions; for de- 
tails, see Ref. [49]. The damaged resistors were replaced 
after Run-1, which significantly reduces the dominant 
contribution to Cpa and the overall magnitude of muon 
losses. 


MAGNETIC FIELD DETERMINATION 


A suite of pulsed-proton NMR probes, each optimized 
for a different function in the analysis chain, measures 
the magnetic field strength [50]. Every ~3 days during 
data taking, a 17-probe NMR trolley [63] measures the 
field at about 9000 locations in azimuth to provide a set 
of 2D field maps. 378 pulsed-NMR probes, located 7.7 cm 
above and below the storage volume, continuously mon- 
itor the field at 72 azimuthal positions, called stations. 
The trolley and fixed probes use petroleum jelly as an 
NMR sample. The probe signals are digitized and ana- 
lyzed [64] to extract a precession frequency proportional 
to the average magnetic field over the NMR sample vol- 
ume. A subset of probes is used to provide feedback to 
the magnet power supply to stabilize the field. 

Calibration procedure fcap: The primary calibra- 
tion uses a probe with a cylindrical water sample. Cor- 
rections are required to relate its frequencies to the pre- 
cession frequency expected from a proton in water at 
the reference temperature 34.7°C. Studies of the cali- 
bration probe in an MRI solenoid precisely determine 


corrections for sample shape, temperature, and magne- 
tization of probe materials to an uncertainty of 15 ppb. 
Cross-calibrations to an absolute ?He magnetometer [65] 
confirm the corrections to better than 38 ppb. 

The calibration probe is installed on a translation stage 
in the SR vacuum. We repeatedly swap the calibration 
probe and a trolley probe into the same location, com- 
pensating for changes of the SR field. This procedure 
determines calibration offsets between individual trolley 
probes and the equivalent w, values. The offsets are due 
primarily to differences in diamagnetic shielding of pro- 
tons in water versus petroleum jelly, sample shape, and 
magnetic perturbations from magnetization of the ma- 
terials used in the probes and trolley body. The trolley 
probe calibration offsets are determined with an average 
uncertainty of 30 ppb. 

Field Tracking (w;,(x,y,@)): The 14 Run-1 trolley 
field measurements bracket muon storage intervals tj, to 
tk+ı- They provide a suite of 2D multipole moments 
(dipole, normal quadrupole, skew quadrupole, ...), which 
the fixed probes track. The fixed probes provide five 
independent moments (four moments for some stations) 
that track the field over 5? in azimuth for each station. 
'The trolley moments are interpolated for times between 
the trolley runs, and the fixed probes continuously track 
changes to five lower-order moments [50]. The fixed 
probe and trolley measurements are synchronized when 
the trolley passes, averaged over each 5? azimuthal seg- 
ment. The trolley run at time t,41 yields a second set of 
moments m'‘"(t,41). The fixed probe moments m?(t, o) 
are used to interpolate the field during muon storage be- 
tween the trolley runs. The uncertainty on the inter- 
polation is estimated from both the k and k + 1 maps 
and a Brownian bridge random walk model. The pro- 
cedure produces interpolated storage volume field maps 
wy (x, y, $) in terms of the equivalent shielded proton fre- 
quency throughout the Run-1 data-taking periods. 

Muon weighting (M(x,y, ¢)): Averaging of the mag- 
netic field weighted by the muon distribution in time 
and space uses the detected positron rates and the muon 
beam distribution measured by the trackers. The inter- 
polated field maps are averaged over periods of roughly 
10s and weighted by the number of detected positrons 
during the same period. The SR guide fields intro- 
duce azimuthal dependencies of the muon distribution 
M(z,y,Q). We determine the muon-weighted average 
magnetic field by summing the field moments m; multi- 
plied by the beam-weighted projections k; for every three- 
hour interval over which the tracker maps and field maps 
are averaged. Along y, the beam is highly symmetric and 
centered, and the skew field moments (derivatives with 
respect to y) are relatively small. The azimuthally aver- 
aged centroid of the beam is displaced radially, leading to 
relative weights for the field dipole, normal quadrupole, 
and normal sextupole of k; — 1.0,0.15, and 0.09, respec- 
tively. An overlay of the azimuthally averaged field con- 


tours on the muon distribution is shown in Fig. 3. The 
combined total uncertainty of a, from probe calibrations, 
field maps, tracker alignment and acceptance, calorime- 
ter acceptance, and beam dynamics modeling is 56 ppb. 


By and By: Two fast transients induced by the dy- 
namics of charging the ESQ system and firing the SR 
kicker magnet slightly influence the actual average field 
seen by the beam compared to its NMR-measured value 
as described above and in Ref. [50]. An eddy current in- 
duced locally in the vacuum chamber structures by the 
kicker system produces a transient magnetic field in the 
storage volume. A Faraday magnetometer installed be- 
tween the kicker plates measured the rotation of polar- 
ized light in a terbium-gallium-garnet crystal from the 
transient field to determine the correction Bz. 


The second transient arises from charging the ESQs, 
where the Lorentz forces induce mechanical vibrations 
in the plates that generate magnetic perturbations. The 
amplitudes and sign of the perturbations vary over the 
two sequences of eight distinct fills that occur in each 
1.4s accelerator supercycle. Customized NMR probes 
measured these transient fields at several positions within 
one ESQ and at the center of each of the other ESQs to 
determine the average field throughout the quadrupole 
volumes. Weighting the temporal behavior of the tran- 
sient fields by the muon decay rate, and correcting for the 
azimuthal fractions of the ring coverage, 8.5% and 43% 
respectively, each transient provides final corrections Dj, 
and B, to a, as listed in Table II. 
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FIG. 3. Azimuthally averaged magnetic field contours 


wp(x, y) overlaid on the time and azimuthally averaged muon 
distribution M (x, y). 


Run wa /2 [Hz] co, /2n [Hz] Rij. x 1000 
la 229081.00(28)  61791871.2(7.1)  3.7073009(45) 
1b 229081.40(24) —61791937.8(7.9) —3.7073024(38) 
le 229081.26(19)  61791845.4(7.7)  3.7073057(31) 
1d 229081.23(16) —61792003.4(6.6) —3.7072957(26) 
Run-i 3.7073003(17) 


TABLE I. Run-1 group measurements of wa, &, and their 
ratios Ri multiplied by 1000. See also Supplemental Mate- 
rial [66]. 


COMPUTING a, AND CONCLUSIONS 


Table I lists the individual measurements of wa and 
Üp» inclusive of all correction terms in Eq. 4, for the four 
run groups, as well as their ratios, R}, (the latter multi- 
plied by 1000). The measurements are largely uncorre- 
lated because the run-group uncertainties are dominated 
by the statistical uncertainty on wa. However, most sys- 
tematic uncertainties for both wa and On measurements, 
and hence for the ratios R}, are fully correlated across 
run groups. The net computed uncertainties (and cor- 
rections) are listed in Table II. The fit of the four run- 
group results has a x?/n.d.f. = 6.8/3, corresponding to 
P(x?) = 7.8%; we consider the P(x?) to be a plausible 
statistical outcome and not indicative of incorrectly esti- 
mated uncertainties. The weighted-average value is Ri, 
= 0.0037073003(16)(6), where the first error is statistical 
and the second is systematic [67]. From Eq. 2, we arrive 
at a determination of the muon anomaly 


a (FNAL) = 116 592 040(54) x 1071! (0.46 ppm), 


where the statistical, systematic, and fundamental con- 
stant uncertainties that are listed in Table II are com- 
bined in quadrature. Our result differs from the SM value 
by 3.30 and agrees with the BNL E821 result. The com- 
bined experimental (Exp) average|[68] is 


a,(Exp) = 116592061(41) x 1071! (0.35 ppm). 


The difference, a, (Exp) — a,,(SM) = (251+ 59) x 10^! , 
has a significance of 4.26. These results are displayed in 
Fig. 4. 

In summary, the findings here confirm the BNL exper- 
imental result and the corresponding experimental aver- 
age increases the significance of the discrepancy between 
the measured and SM predicted a, to 4.20. This result 
will further motivate the development of SM extensions, 
including those having new couplings to leptons. 

Following the Run-1 measurements, improvements to 
the temperature in the experimental hall have led to 
greater magnetic field and detector gain stability. An 
upgrade to the kicker enables the incoming beam to be 
stored in the center of the storage aperture, thus reducing 
various beam dynamics effects. These changes, amongst 
others, will lead to higher precision in future publications. 


Quantity Correction terms Uncertainty 

(ppb) (ppb) 
wa (statistical) = 434 
wa’ (systematic) - 56 
Ce 489 53 
Cp 180 13 
Cmi -11 5 
Coa -158 75 
fealib (wp (x, y, à) x M(x,y, $)) 7 56 
By, -27 37 
By -17 92 
My (34.7°)/pHe = 10 
m, /me - 22 
Ge/2 - 0 
'Total systematic E 157 
Total fundamental factors = 25 
Totals 544 462 


TABLE II. Values and uncertainties of the Ri, correction 
terms in Eq. 4, and uncertainties due to the constants in Eq. 2 
for ap. Positive C; increase a, and positive B; decrease ap. 
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FIG. 4. From top to bottom: experimental values of a, 
from BNL E821, this measurement, and the combined aver- 
age. The inner tick marks indicate the statistical contribution 
to the total uncertainties. The Muon g — 2 Theory Initiative 
recommended value [13] for the standard model is also shown. 
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